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The onset of the PF2 superfluid phase in high-density neutron matter is studied within the BCS 
framework with two and three body forces. Owing to the strong correlations the energy gap is so 
sizeably quenched as to demand to reconsider the role of superfluidity in neutron-star phenomena. 
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I. INTRODUCTION 

So far, the neutron stars (NS) have been considered as 
a rich laboratory of various superfluid phases of nuclear 
matter p]-[lOj. Recently the interest has been focussed 
on the NS interior, where both the vortex pinning re- 
sponsible for the observed period glitches [llj and the nu- 
cleon superfluidity responsible for the main cooling mech- 
anisms |12l | are supposed to be located. In particular the 
recent observations of cooling in the NS of Cassiopeia A 
have been considered to be a direct evidence of the ^PF2 
pairing in NS core, and the energy gap needed to explain 
the data has been estimated to be around 0.1 MeVT^I- 
From the theoretical viewpoint, the onset of the ^PF2 
superfluid state occurs in regions above the saturation 
density. This event rises the problem of the pairing for- 
mation in high density nuclear matter, that however all 
previous calculations neglect (see Ref.[13| and references 
therein). 

The density region, where the '^PF2 component of the 
nuclear interaction is the most attractive, may extend 
to several times the saturation density of nuclear mat- 
ter. In such a high-density regime, the short-range cor- 
relations are so strong that the momentum distribution 
around the Fermi level departs significantly from the typ- 
ical profile of a degenerate Fermi system. This departure 
is measured by the so-called Z- factor (0 < Z < 1) [l4| . 
Since the deformation of the Fermi surface hinders par- 
ticle transitions around the Fermi energy ep, the pair- 
ing gap is expected to get reduced in any case, even for 
the low density ^5*0 pairing channel (see, e.g. ,[15]). Con- 
cerning the high-density neutron-neutron pairing in the 
^PF2 channel, the calculations so far reported in the lit- 
erature ignore the effect of the Z-factor, owing to the 
large uncertainty still existing in the magnitude of the 
nucleon-nucleon (TV TV) interaction in the ^PF2 channel, 
including both the two body force (2BF) and the three- 
body force (3BF) component. Other many-body effects, 
such as screening effect, have been neglected for the same 
reason. 

In this note, we present a study of the Z-factor effect 
on the ^PF2 pairing in pure neutron matter. In principle, 
we should consider asymmetric nuclear matter for appli- 
cation to NS core, but the small proton fraction is not 
relevant in this context, as discussed below. The defor- 



mation of the Fermi surface and the Z-factor are studied 
in the framework of the Brueckner theory with 2BF and 
3BF [i^. The energy gap is then calculated within the 
generalized BCS theory [l7], including in the pairing in- 
teraction not only 2BF, but also 3BF. The latter, in fact, 
is dominant at high density and therefore it is expected to 
directly influence pairing gap in addition to the Z-factor. 



II. FORMALISM AND RESULTS 
A. Nucleon propagator in neutron matter 

The neutron Green's function is given by 
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where ep denotes the Fermi energy and S(p,(u;) the self- 
energy. Expanding the latter in a series of powers of the 
quasiparticle energy around the Fermi surface, we obtain 
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where the quasi-particle energy and the quasi-particle 
strength are 
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respectively. The quasiparticle strength Z(p) measures 
the deviation of a correlated Fermi system from the ideal 
degenerate Fermi gas. The occupation numbers n(p) and 
the Z{p) factors have been calculated in the framework 
of the Brueckner theory [l^ [13] . Grace to the inclusion 
of the three body forces the hole-line expansion can be 
extended up to high densities. The self-energy, truncated 
to the second order (see Ref. '2&\ for details), provides us 
with a g ood reproduction of the empirical nuclear mean 
field [2l| and the optical-model potential [22| , so that we 
are quite confident that the next orders are irrelevant for 
the present calculation. We employ as 2BF the meson 
exchange Bonn B potential [23[, whose meson param- 
eters are constrained by the fit with the experimental 



phase shifts of NN scattering. The microscopic meson 
exchange 3BF is the one constructed by Li et al. p^ . 
It is consistent with the 2BF because it adopts the same 
meson parameters as Bonn B, so that there are no free 
parameters in the model. 

The Z{p) factor is related to the depletion of the oc- 
cupation number n{jp) around the Fermi surface. Ac- 
cording to the Migdal-Luttinger theorem [24| its value 
Z — Z(p^)(Z-factor) equals the discontinuity of the oc- 
cupation number at the Fermi surface, i.e., 
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where p^ is the Fermi momentum. In our approximation 
Yj{p,ijj) = I]i(p, oj) + I]2(p, tj), where I]i(p,aj) determines 
the left limit and Yi2{p,ijj) the right limit of the preceding 
equation for e — ;► 0. 



the ideal Fermi gas (solid line) due to the strong short- 
range correlations. As expected, the deviation is slightly 
enhanced by 3BF. In Fig. 1(b) the calculated Z-factor is 
displayed vs. density in the two different approximations 
for the self-energy, i.e., E « Ei and E « Ei + E2, re- 
spectively. The calculation of Zp from Eq.(3) requires 
a high numerical accuracy: increasing the accuracy the 
calculated Zp gets lower and lower until the converging 
value is reached. It is noticed that without 3BF, the Z- 
factors decrease slowly as a function of density. Adding 
the contribution of E2 leads to an overall reduction of the 
Z factor. The 3BF reduces further the Z-factor and its 
effect increases rapidly with density. As a consequence, 
including 3BF makes the decrease of the Z-factor as a 
function of density much rapid than that obtained by 
adopting pure 2BF. Therefore 3BF induces a strong ex- 
tra deviation from the ideal Fermi gas model. 
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B. Gap equation in the ^PF2 channel 

The ^PF2 superfluidity in pure neutron matter has 
been investigated by using various theoretical approaches 
|25l - l3l| with 2BF, and later extended to microscopic 3BF 
forces by Zuo et al. [13] ■ In this case, the pairing gaps 
are determined by the two coupled equations: 
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FIG. 1: (Color Online) Occupation numbers vs. momentum 
at a given density (upper panel) and Z-factors vs. density 
(lower panel) in pure neutron matter. The effect of 3BF is 
shown in both panels. In the lower panel the calculations are 
reported for two approximations of the self-energy 

In Fig. 1(a), we display the calculated occupation 
probability in pure neutron matter at density p = 0.3 
fm^'^. One easily observes the remarkable deviation from 



where E^ = (e^ - /z)^ + A^ and A^ 
Vl,l'{p,p') are the matrix elements of the realistic A^A^ 
interaction in the coupled ^PF2 channel. In the gap equa- 
tion, the Z-factors and the single-particle energy tp are 
calculated from the Brueckner theory. As for the pair- 
ing interaction Vp.L'iPiP') in the ^PF2 coupled channel, 
we adopt the same 2BF and 3BF as in the Brueckner 
calculation. The 3BF cannot be neglected in the gap cal- 
culation because the ^PF2 pairing is expected to occur 
in the high density domain, where 3BF is quite sizeable 
specially in the ^PF2 channel. 

The results are summarized in Fig. 2. Neglecting the 
Z-factor effect (upper panel), the magnitude of the ^PF2 
gap with the new interaction dos not differ from previ- 
ous calculations with AV18 potential 17i]. On the other 
hand, the 3BF enhances the ^PF2 superfluidity signif- 
icantly at higher densities. But, as shown in the lower 
panel of Fig. 2, the introduction of the Z-factor effect fa- 
tally quenches the pairing gaps to a value less 0.05 MeV, 
one order of magnitude smaller than the value with full 
interaction. The effect of the Z-factor appears to be 
extremely sizeable at high densities. It is worth notic- 
ing that this effect is opposite to the 3BF effect on the 
^PF2 pairing in neutron matter: the former turns out 
to be much stronger than the latter. In conclusion, the 
departure of the system from the pure degenerate limit 



drives the pairing attenuation at high density, regardless 
of whether the 3BF is included or not. 
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III. CONCLUSIONS 

In conclusion, we have studied the anisotropic ^PF2 
pairing in pure neutron matter. The effects of the Fermi 
surface depletion (Z-factor) have been included in the 
calculation of the energy gap. In the pure degenerate 
limit, the ^PF2 superfluid phase extends over a broad 
density range with a gap peak value of about 0.2 MeV 
without 3BF and 0.5 MeV with 3BF. The inclusion of the 
Z-factor leads to a rapid decrease of the gap magnitude 
by one order of magnitude: its peak value drops to less 
than 0.05 MeV and the superfluidity domain shrinks to 
0.1 — 0.4 fm~^. In neutron stars the proton fraction in 
/3-equilibrium with neutrons could in principle affect the 
'^Pi^2 pairing, but in the non- vanishing gap range it is less 
than 15% of the total density. Recent calculations on the 
isospin dependence of the quasi-particle strength show 
that the enhancement of the neutron Z-factor is negligi- 
ble for such small proton fraction 13311. Even including 
the additional screening suppression [l^, the conclusion 
remains that the departure from the Fermi gas limit is 
the main cause of pairing disappearance in high-density 
nuclear matter. This result makes doubtful the role of 
the ^PF2 pairing in NS core. 



FIG. 2: (Color Online)Effect of the Fermi surface depletion on 
''PF2 pairing gap in pure neutron matter. Notice the y-scale 
change from (a) to (b) 



In Ref . [3^ , the effect of the Z-factor was investigated 
also for the ^5*0 pairing in neutron matter, and it was 
shown to reduce the energy gap. Comparing with the 
results of Ref.[32'|, the same effect on the '^PF2 superflu- 
idity appears to be much stronger. This is simply due to 
the fact that in the density region of the '^PF2 superfluid- 
ity the deviation of neutron matter from the degenerate 
Fermi gas limit becomes much larger. 
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